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Abstract. Anomalous transport processes in which the 
variance of the distance travelled does not necessarily in- 
crease linearly with time: (Acc^) cx with < a < 2 are 
modelled using the formalism of continuous time random 
walks. We compute particle propagators which have the 
required dependence on space and time and use these to 
find the spatial dependence of the synchrotron radiation 
emitted by a population of continuously injected electrons. 
As the electrons are transported away from the source 
they cool, and the synchrotron spectrum softens. Sub- 
diffusive (a < 1) transport - corresponding to stochastic 
trapping, or restriction of the transport across the aver- 
age direction of a stochastic magnetic field - produces a 
much slower rate of change of spectral index than does 
supra-diffusion (1 < a < 2) - which occurs when particles 
move almost without scattering, in a field containing large 
ordered regions. Application to the diffuse emission of the 
outer parts of the Coma cluster favours an interpretation 
in terms of supra-diffusion. 

Key words: diffusion - magnetic fields - plasmas - turbu- 
lence - galaxies: clusters: Coma; radio continuum: galaxies 



1. Introduction 

Many of the properties of cosmic rays can be understood 
on the basis of a propagation model in which they exe- 
cute isotropic diffusion between injection at sources dis- 
tributed in the galactic disk and escape at certain bound- 
aries whilst possibly being carried out of the galaxy in a 
wind (Lerche & Schlickeiser 1982 Bcrezinskii et al. 1990| , 



Webber et al. 1992, Blocmcn ct al. 1993). On the other 



hand, it has been known for many years that these parti- 
cles are 'magnetised' in the sense that their Larmor radius 
is much smaller than the mean free path between 'col- 
lisions' implied by the diffusion coefficient. This means 
that particles diffuse primarily along the magnetic field. 



and only very slowly across it. Since the galactic magnetic 
field lies predominantly in the galactic plane (e.g., Zweibel 
& Heiles 1997), it is surprising that isotropic diffusion de- 



scribes cosmic ray propagation so well. One might instead 
expect cosmic rays to be transported anomalously slowly 



in the direction normal to the disk (Getmantsev 1963). 
The success of the simple diffusion model of propaga- 
tion (or even of the more primitive but related 'leaky- 
box' model) appears to show that this is not the case, at 
least for the bulk of cosmic rays. This is presumably due 
to fluctuations in the magnetic field, possibly connected 
with instabilities, which cause a given field line to wan- 
der out of the galactic plane, and give rise to an effective 
diffusion across t he dir ection o f the mean magnetic field 
(Jokipii & Parker |l969| , Jokipii |l973|) . 

Nevertheless, it may not always be the case that the 
transport process can be described by an effective diffusion 
coefficient. Since its first mention by Getmantsev (1963), 
anomalous transport, i.e., transport in which the mean 
square deviation of a particle from its position at a given 
time (Ax^) does not depend linearly on the elapsed time, 
has been discussed in two astrophysical contexts. Firstly, 
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Chuvilgin & Ptuskin (1993) derived a kinetic equation for 
cosmic ray transport for the case in which (Ax^) cx t^^"^, 
but pointed out that this type of transport can be ex- 
pected only when considering timescales short compared 
to the time taken for a particle to effectively decorrelate 
from a given magnetic field line. Many effects such as par- 
ticle drifts, temporal changes in the magnetic field or just 
the chaotic structure of the field itself can cause decor- 
relation, but although it is very difhcult to estimate this 
time in an astrophysical situation, the success of the dif- 
fusion model suggests that cosmic rays do indeed decorre- 
late from the field faster than they escape from the galaxy 
(which takes about 10* years). 

The second astrophysical context is that of particle 
acceleration at a shock front. In this case, the natural 
timescale is the acceleration time, which is thought to be a 
strong function of energy, and to vary over several orders of 
magnitude for those cosmic rays accelerated at supernova 
shocks (Dendy et al |l995i Duffy et al. [19951) . Here too, it 



2 



B. R. Ragot et al.: Anomalous transport of cosmic ray electrons 



is difficult to make a realistic estimate of the decorrela- 
tion time of a particle from the magnetic field. However, 
predictions can be made of the spectrum and spatial dis- 
tribution of those particles undergoing anomalous trans- 
port, which might enable them to be distinguished from 



diffusing particles (Kirk et al. 1996a) 



A problem which is similar in many respects arises 
in the confinement of plasma in fusion devices, and in 
this context there has been considerable interest in recent 
years in anomalous, non-diffusive transport models (Roch- 
ester & Roscnbluth|l978| Kadomtsev fc Po gutse |1979| , 
IsichenkopMla|,[l991b|, Rax & White |1992|, Wang et 



al. 1995). In particular, use of the formalism developed for 
continuous time random walks (CTRW - see MontroU & 



Weiss 1965) has been advanced (Balescu 1995). This ap- 
proach is widely used in the problem of anomalous trans- 
port in random media (e.g., Bouchaud & Georges 1990] , 
Shlesinger et al. 19931 ). In this paper, we discuss these 
techniques and show how they relate to the work already 
done on anomalous transport in astrophysics. We point 
out that they offer a more general approach to the prob- 
lem, enabling one to relax the rather restrictive assump- 
tions concerning the statistical properties of the magnetic 
field employed hitherto. As an example, we examine a 
new astrophysical application - that of the synchrotron 
emission of a population of relativistic electrons which 
lose energy whilst being transported away from the site 
of their acceleration. Both the spatial extent of the ra- 
diating electrons and the observed spectrum depend on 
the nature of the transport process. These calculations 
are especially relevant to the interpretation of high res- 
olution radio observations in several astrophysical situa- 
tions, including spiral galaxies seen 'edge-on' (e.g.. Hum- 
mel et al. 1991) as well the diffuse emission from clusters 
of galaxies (Schlickeiser et al. 1987, Kirk et al. 1996b). 



2. Anomalous transport 

In a regular magnetic field topology, the transport of 
charged particles across the field is due to the collisions 
of the particles and their finite Larmor radii. However, a 
perturbation of the magnetic field results in a wandering 
of the field lines and a potentially much faster transport 
of the particles. This transport is widely attributed to the 
effect of the large-scale random field component, which 
would produce, following the quasilinear theory, a diffu- 
sion of the field lines across the direction of the the average 



field (Jokipii & Parker 1969, Kadomtsev & Pogutse 1979) 



In fusion applications, it is generally assumed that particle 
collisions lead to diffusive transport along each individual 
field line. In astrophysical applications, where collisions 
can be neglected, it is usual to assume that small scale 
fluctuations in the magnetic field play this role, so that 
here too , par ticles diffuse along field lines (Chuvilgin & 
Ptuskin 1993). As long as the particles remain correlated 
to a given patch of field lines, the combination of these 



two diffusions results in sub-diffusion of the particles i.e., 
{Ax'^{t)) oc t", with a = 1/2, as described by Getmant- 
sev ( 1963| ). (In the general case, sub-diffusion refers to all 
such processes when < a < 1.) However, the exponen- 
tial divergence of neighbouring field lines and the resulting 
stretching of a field line patch lead to decorrelation of a 
particle from its field line, and thus to the large-scale dif- 
fusion of the part icles, as pointed out by Rechester and 
Rosenbluth ( 1978|) , with a diffusion cocfhcient given by 
the expression: 



(1) 



where IcS = log(l/fco'5), 5 — lc\/ k^/ where Ic is the 
exponentiation length of the field lines. Here Dst is the 
quasilinear diffusion coefhcicnt of the field lines, which 
has the dimensions of a length, is a characteristic wave 
number of the magnetic turbulence, Ky and Kj^ are the 
quasilinear diffusion coefficients of the particles along and 
across the direction of the local magnetic field, respec- 
tively. The transport of particles across a turbulent mag- 
netic field is thus sub-diffusive for short timescales, but 
crosses over to normal diffusion in the limit of the long 
timescales. In the following, we investigate the anomalous 
transport regime, which applies if the natural timescale 
of a particular problem - such as escape from the galaxy, 
acceleration at a shock, or loss of energy by synchrotron 
radiation, is shorter than the decorrelation time. 

In a statistically homogeneous medium, the density of 
particles, n{x,t), is related to the source Q{x',t') by the 
propagator P: 



(2) 



For diffusive transport, P typically contains a factor 
exp(— |a;p/_Di), and is the Green's function of the diffusion 
equation. 

Fo r sub-diffusive transport with a = 1/2, Rax & White 
(1992) have determined the propagator by combining two 
diffusive propagators using a Wiener integration method. 
They considered the transport in cylindrical symmetry 
across the z direction (i.e., a 2-dimcnsional problem, with 
\z\ playing the role of time), and obtained: 



P(r, t) 



Hit) 



27r {Ar^Dl,K„t) 



1/3 



exp 



, 3^4/3 



1/3 



(3) 



where r is the radius in cylindrical coordinates, and H{t) 
the Heaviside function. Duffy et al. (1995) considered the 
transport perpendicular to a plane shock front. In a carte- 
sian coordinate system with the shock in the y-z plane, 
and assuming no gradients in the y direction are present, 
they obtained the propagator 

^ f exp [-xy (iPMlzl) - zV (4aC|,0] 



(4) 



B. R. Ragot et al.: Anomalous transport of cosmic ray electrons 



3 



Approximating the integral using the method of steepest 
descents, they found 



p{x,t) « v{DMKy'\x\f^T'^' 



exp 



-/3|a;|4/3 



r)2/3 1/3,1/3 



(5) 



with T] = 2^/3 (Stt) and /? = 3/2^/3^ For this 1- 
dimensional prob lem of a = 1/2 sub-diffusion, Chuvilgin 
& Ptuskin ( 1993 ) derived an equation describing the evo- 
lution of the particle density (their Eq. B.12) and solved it 
to find the above propagator. E ssenti ally the same equa- 
tion was also found by Balescu ( 1995 ). It can be written: 



dtn{x,t) = Do 



dT 



l/rr 



T 



3/2 



V^n{x, t — t) 



(6) 



with Dq — a'^/2TD- This is a non-Markovian diffusion 
equation: the integral term is characteristic of the long- 
time memory of the dynamics. 

In a realistic situation, the topology of the magnetic 
field may be more complicated than just a pure stochastic 
sea with its associated diffusion of the field lines. In fusion 
plasmas, for example, there can exist ordered structures, 
'stability islands', in which particles can be trapped for 
long periods of time, (referred to as 'sticking'), leading to 
sub-diffusive large-scale transport of the field lines them- 
selves (e.g., White et al. 1993D . On the other hand, the 
field lines may also wander faster than implied by dif- 
fusion. In this case, the field lines are said to perform 
'flights', during which they maintain almost the same di- 
rection for a relatively long distance. Ultimately, the large- 
scale transport of field lines is the result of competition 
between 'sticking' and 'flights', and can yield transport 
regimes of many different kinds, ranging from slow sub- 
diffusion (almost perfect sticking) to fast supra-diffusion 
(domination of flights). In terms of a, the physically rel- 
evant range is < a < 2, with a = 2 corresponding to 
transport completely dominated by a single straight-line 
flight. 

To find the particle propagator, it is necessary to com- 
bine the propagator for field lines with that for particle 
motion along the field, as described above for the case of 
a = 5. Here too, the assumption of diffusive transport can 
be generalised. Particles may, in fact, undergo no scatter- 
ing at all, in which case they propagate ballistically along 



the field lines, as assumed by Achterberg & Ball (1994) 



On the other hand, they may be trapped between mag- 
netic mirrors on a segment of a field line. Once again, the 
transport can be characterised by an index a which lies 
between and 2. However, as we show below, the com- 
bination of two propagators does not extend the overall 
range of a which is permitted. In fact, the effect of su- 
perimposing transport with the two indices ai and a2 is 
described by a single value a — aiQ;2/2. 



3. Continuous time random walks 

We do not aim to describe the microscopic dynamics of 
particles moving in a turbulent magnetic field in a realistic 
way, but rather look for a model which is able to repro- 
duce their transport characteristics and gives the density 
profile. In terms of the propagator P{x,t), the density is 
given by Eq. (2), once the rate Q{x,t) at which particles 
are injected is specified. This can then be used to compute, 
for example, the synchrotron radiation from electrons in- 
jected at a plane surface (such as a spiral galaxy) and 
subsequently transported outwards. 

The method we adopt to model P is that of contin- 
uous time random walks (CTRW). In this, the transport 
of a particle is modelled as a succession of instantaneous 
jumps of arbitrary length, separated by pauses of arbitrary 
duration. Each of these steps is independent of the previ- 
ous steps. We denote by f{x) the probabihty distribution 
function (PDF) of a jump described by x, and by ip{t) the 
PDF of a pause of duration t, and use their Fourier and 
Laplace transforms: 



/>)= / dxf{x)e' 



+ 00 



dt'ip{t)e- 



(7) 



The type of anomalous transport characterised by a 
and discussed above arises from the asymptotic behaviour 
of P{x,t) as \x\ — + +00 and t +00. This corresponds 

to the |fc| ^ 0, s ^ asymptotic behaviour of P{k,s), 
for which the expansions of -0(5) and /(fc) around are 
required. To keep the treatment general, we consider first 
the motion of a particle in a space of dimension d. The 
functions i/'(s) and f{k) are then prescribed as 



V3(s) 
/(fc) 



TjjS 





Ifel 



0, 



(8) 
(9) 



with < a, /3 < 2, and td, a positive constants. In Eq. (9), 
we have assumed that / depends only on |A;|, i.e., that 
there is no preferred direction for a particle jump. The 
particular values a = 1, /3 = 2 yield a diffusive process 
with a Gaussian propagator, whereas the case /3 < 2, leads 
to a random walk with an infinite second moment (x^), 
corresponding to Levy flights. 

In a classical paper, MontroU and Weiss ( |1965 ) have 
shown that the Fourier-Laplace transform of P{x,t) has 
the following form 



P{k,s) 



1 - ij{s) 



.[l-0(5)/(fc)] 



(10) 



Inserting the expansions Eq. (9) we find that for large \x\ 
and t 



dx \x\'^P{x,t) 
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dfc 



Ax 



7(2d) 



(11) 



and the transport coefficient is n = 2a/ (3. In the follow- 
ing, we find it convenient for the inversion of the Fourier 
transforms to choose (3 — 2, and describe all transport 
regimes using the parameter range < a < 2: 



il){s) = 1 - T^s" , < a < 2 



(12) 
(13) 



for small s and |fc|. The choice a = \ has the same asymp- 
totic behaviour as the function ^ = cxp(— s/t/j), i.e., 
= S{t—TD), which corresponds to scattering after a fixed 
time interval, yielding diffusion. However, for a ^ 1,2, we 
have 



1 



t 



for t — !• +00. This is a long-tailed distribution function. 
Its slow decrease at large t leads to the same kind of mem- 
ory effect built in to the non-Markovian transport equa- 
tion (6) constructed by Chuvilgin & Ptuskin ( 1993 ) and 
Balescu (|l995|) for the case a < 1. We do not consider 
a > 2, since in this case the quantity {x^) /t^ increases 
without limit for large t, and would at some stage exceed 
the actual physical particle speed. 

From the expansions (13) the asymptotic behaviour 
of the propagator at large |a;| and t can be found by 
an inverse Fourier transformation followed by an in- 
verse Laplace transformation, which is performed approx- 
imately using the method of steepest descents. The re- 
sulting expression can be used for all values of x and t 
once the normalisation has been corrected. It represents 
a convenient one-parameter model of particle transport 
which exhibits the anomalous behaviour in which we are 
interested. The details of the derivation are given in Ap- 
pendix A. Here we simply present the results, expressed 
in terms of a dimensionless similarity variable defined ac- 
cording to 



x\ /rD\"/2 



Using this variable, the propagator separates: 
and we find 

S,(0 = C,r'''-"^/^'-"^exp[-A^2/(2-«) 



(15) 



(16) 



(17) 



where the (positive) constants Cd and A are given 
in the appendix (Eqs. A3 and A4). These propagators 
are normalised such that particle number is conserved: 
/ d'''xP{x,t) — 1. They display at all t the behaviour re- 
quired of the asymptotic dependence of the variance in 
anomalous transport. In particular, they possess for all t 
the property 



(18) 



For one-dimensional propagation (d — 1), which de- 
scribes, for example, the transport of particles in x away 
from a uniform source in the y~z plane, we have for the 
standard diffusive case {a — 1) the well-known form 



Pix,t) 



1 



exp [-eV4] 



(19) 



with D = a'^/{2TD) and ^ = In the a = 1/2 

sub-diffusive case, which can be considered as a double 
diffusion, we find: 



(14) Pix,t) 



1 



t 



1/4 



-1/3 



exp 



^3(e/4) 



4/3 



(20) 



which agrees with the propagator for d = 1 given in 



Eq. (5) and derived by Duffy et al. ( |1995D , with a = 
{2DM?'^k\'^/u^'^ and td = (2Dm)'/^/^^^/'. (Sim- 
ilarly, the case d = 2 reproduces the propagator (3) of 
Rax & White 1992). In analogy with the formulation of 



Eq. (4), it is possible to consider other values of a as aris- 
ing from a convolution of a propagator describing motion 
along the field line together with one describing the wan- 
dering of the field. Using the method of steepest descents, 
it is straightforward to show that Pd,a has the following 
convolution property: 



+ OC 



dsPd,^i(a;, |s|)Pi,P2(s,i) = Pd,A.iM2/2(a^:0 , (21) 



(note that the normalisation is preserved) so that the com- 
bined transport process can always be described using a 
value of a between and 2. 

To illustrate the properties of these propagators, they 
are plotted for three different values of a in Fig. (|l|). The 
diffusive propagator is simply a gaussian curve given by 
Eq. (19) according to which the probability of finding a 
particle at a particular distance drops off smoothly with 
increasing distance from the point (plane) of injection. 
The sub-diffusive propagator, on the other hand, falls off 
much more steeply initially (it has an integrable singular- 
ity at the point of injection), but at larger distances de- 
creases more slowly than diffusion. Acting on an ensemble 
of particles, sub-diffusive transport tends to confine some 
of them close to the injection point, but propagates others 
very rapidly to large distance. This is a result of the larger 
spread in position of sub-diffusing particles of a given age. 
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Fig. 1. The dimensionless propagator (Eq. 16) for 
one-dimensional transport in three different regimes: standard 
diffusion {a — 1, solid line), sub-diffusion (a = 1/2 dotted 
line) and supra-diffusion (a = 3/2 dashed line), as a function 
of the dimensionless similarity variable ^, which, for fixed time, 
is proportional to distance (see the definition in Eq. 15). 



compared with diffusing ones. Thus, despite the fact that 
sub-difFusion produces, on average, slower transport, a mi- 
nority of particles experiences comparatively rapid trans- 
port. Supra-difFusion displays the opposite behaviour and 
is similar to ballistic transport, or propagation at con- 
stant speed without scattering. The propagator is strongly 
peaked around the value ^ « 1. Very few particles are to 
be found lingering close to the point of injection, and very 
few escape to large distance. 

Equation (17) is the basic result of this section. We 
now turn to the question of the observable consequences 
of these propagators. 



4. Synchrotron radiation from anomalously trans- 
ported electrons 

The synchrotron radiation emitted by relativistic electrons 
is the most important diagnostic available for the study of 
the transport of these particles in astrophysical environ- 
ments. In this section we compute the spatial dependence 
of the surface brightness and of the spectral index which is 
to be expected when the different kinds of particle trans- 
port discussed above dominate. 

The energy losses suffered by a relativistic parti- 
cle emitting synchrotron radiation, or undergoing inverse 
compton scattering whilst propagating through a homo- 
geneous magnetic field, or a homogeneous soft radiation 



field can be described by the equation 
d7 2 

dT = 



(22) 



where 7 is the Lorentz factor of the particle and g a factor 
which depends on the magnetic field strength, and the 
energy density of the photon field. In units suited to the 
application to clusters of galaxies (see below) one has 



= 4.1 X 10-14 {Bl^ 



12.3) 



yr 



(23) 



where the photon field is assumed to be that of the cosmic 
microwave background, which is for this purpose equiva- 
lent to a magnetic field strength of 3.5 /iG. Defining the 
differential rate dQ at which electrons are injected into the 
system at time t as dQ = Q{x,j,t)d'^xd'y, the resulting 
electron density, allowing for energy losses, is 



nOO n nt 

/ di d'^x' dt'P{x-x',t-t') 

J ■y J J —OQ 

Q{x',j',t')S (j- 



1 + gi{t - t') 



(24) 



Consider electrons continuously injected with a power-law 
distribution in 7 at the point a; = 0, i.e.. 



Q{x,^,t) - Qo-^-''5[x) 



(25) 



where Qo is a constant. At all points and at all energies of 
interest the distribution has had sufficient time to reach a 
stationary state, which is given by the integral 



n(x,-i) ^ Qo-f- 



1/(97) 



dt'P(x,t'){l-gjt'f-^ . (26) 



In order to compute the surface brightness on the sky of 
the synchrotron radiation from these electrons, two more 
integrations are required - firstly over the line of sight 
through the source and secondly over 7, after multiply- 
ing n by the synchrotron kernel. Furthermore, informa- 
tion about the magnitude and direction of the magnetic 
field along the line of sight is needed. For our purpose, it is 
sufficient to use the delta-function or 'monochromatic' ap- 
proximation to the synchrotron kernel (see, for example, 
Mastichiadis & Kirk 1995| ) and we will simplify the dis- 
cussion by using a constant, angle-averaged emissivity in 
the line of sight integral (see, however, Crusius & Schlick- 
As a result, the surface brightness at a given 



198S 



eiser 

frequency is simply proportional to Jdzrt(a;,7), where z 
is measured along the line of sight and where 7 is fixed, 
being proportional to the square root of the observing fre- 
quency. 

In the one-dimensional case in which electrons are 
injected in a plane, and the line of sight is parallel to 
this plane, displaced by a distance x, the resulting sur- 
face brightness is directly proportional to n(x,7) given 
by Eq. (26). This expression is simple to evaluate numer- 
ically. For various special cases, such as diffusion {a ~ 1), 
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or injection with p = 2 the integral can be written in 
closed form (see Appendix B for details). The resulting 
profile is shown in Fig. |^ for the three qualitatively differ- 
ent types of transport considered in Sect. ^. It is important 




Fig. 2. The surface brightness of synchrotron radiation from a 
planar source of relativistic electrons, seen edge-on, as a func- 
tion of distance x from the plane of injection, normalised to 
unity at x = 0.5. Distance is expressed in units of the scale 
length Xs (Eq. 27). Electrons are injected at a; = with a 
spectrum Q oc 7~^'^ 



to note that the shapes of these profiles are characteris- 
tic of the type of transport involved and are independent 
of the frequency of observation. The assumptions entering 
the calculation are that the magnetic field strength is con- 
stant, that the source has the stated geometry, and that 
the electron injection has been steady over a sufficiently 
long period of time. The abscissa in Fig. || is the scaled 
radius x ~ x/xs, where 



V2(57rD)"/2 



(27) 



which is in general a function of observing frequency. Sub- 
diffusive transport {a < 1), which tends to trap particles 
close to the point of injection, displays an integrable singu- 
larity at cc = 0, whereas supra-diffusive transport, which 
quickly moves particles away from the source, initially 
increases at small x. The brightness profile for diffusive 
transport is monotonically decreasing with X. In Fig. H 
the surface brightness is arbitrarily normalised to unity 
aX X — 0.5. In a realistic situation, in which the electron 
source is distributed in space, it may not be possible to de- 
tect these features. In fact, the smooth almost exponential 
decrease of surface brightness with distance means that it 



would be difficult if not impossible to distinguish obser- 
vationally between the different types of transport from 
these profiles alone. 

A more promising approach lies in the measurement of 
synchrotron spectra as a function of position. In the delta- 
function approximation, the synchrotron spectral index 
as is simply related to the 7-dependence of the electron 
density integrated along the line of sight. To compute this, 
one must take account of the dependence of the scaling 
radius on 7. Defining the frequency dependence of the 
surface brightness by Ii, oc v"" , we have 



2 



1 



d\og{N) 
d\og{x) 



(28) 



where N — J dzn{x,"f). The resulting spectral index as a 
function of x is plotted in Fig. for p = 2.5. 
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Fig. 3. The spectral index of synchrotron emission as a 
function of (scaled) radius. The dotted line corresponds to 
sub-diffusion with a = 0.5, the solid line to diffusion {a = 1) 
and the dashed line to supra-diffusion {a — 1.5). Electrons 
are injected continuously in the plane x = with a spectrum 
Q oc 7~^'^. The total emission integrated over x is thus a pure 
power-law of index = —1.25. 



This figure shows that the different types of transport 
produce very different spectra as the particles propagate 
away from the source. Close to the plane of injection, 
supra-diffusion gives approximately the "uncooled" value 
of the spectral index as = —0.75. This is because only 
recently injected particles remain close to x = - older 
particles have a very small chance of return in the supra- 
diffusive case. As a decreases, more and more older parti- 
cles are confined close to the point of injection, so that the 
spectrum for diffusive and sub-diffusive transport is softer. 
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In general, the age distribution of particles produced by 
each transport process determines the spectrum. Thus, 
supra-difFusion (a = 1.5, dashed line) has a very narrow 
spread of particle age at a given x. Once the distance ex- 
ceeds Xs, the spectrum becomes very soft, since the num- 
ber of radiating particles decreases rapidly. In the case of 
diffusion (q = 1, solid line) the age distribution is gaussian 
and produces a more gradual softening, whereas the sub- 
diffusive case {a = 0.5, dotted line) not only traps some 
particles close to x = but, also transports a few par- 
ticles out rapidly outwards. Thus, the spectrum remains 
rather hard even at large radius. An observationally rele- 
vant measure of this dependence is the change in spectral 
index over the distance required for the brightness to drop 
by one order of magnitude from its value at x = 0.5. For 
supra-diffusion, the spectrum softens by 1.8 from —1 at 
X = —.5 to —2.9 at X — 1.6; for diffusion it softens by 0.8 
from —1.1 at a; = 0.5 to —1.9 at x = 2 and for sub-diffusion 
we find a change of only 0.3 from —1.2 at a; = 0.5 to —1.5 
at a; = 2.2. Thus, a relatively slow change in synchrotron 
spectral index indicates sub-diffusive behaviour. 

As an example, we consider the diffuse radio emission 
observed from the Coma cluster of galaxies. There is ac- 
tive debate as to the origin of the relativistic electrons 
responsible for this emission. It appears that simple diffu- 
sion from one of the galaxies near the centre is not, by it- 
self, a viable explanation (e.g.. Kirk et al. |1996b ) and that 
some kind of distributed inj ectio n and/or acceleration is 
necessary (Schlickeiser et al . |1987 ). However, observations 
by Giovannini et al. ( |1993|) have shown that there is a 
distinction between the central parts of the cluster (say 
within a radius of 500 kpc) and the outer parts (in this 
paper we adopt the value Hq — 50kms~^ Mpc~^, so that 
1 sec of arc is equivalent to 700 pc). If we assume that 
the processes of injection and acceleration occur within a 
sphere of radius 500 kpc, we can apply the above propa- 
gators to model the surface brightness and spectral index 
of the emission in the outer parts. 

Figure ^ shows the observed emission at 430 MHz, 
as presented by Valtaoja ( |1984D , from the data of 
Hanisch ( 198C| ). The horizontal "error bars" in this repre- 
sentation are an expression of spread in size of the emission 
when measured along the E-W and N-S directions. Super- 
imposed on the data are three models computed using the 
one-dimensional propagators. Clearly, the (almost) spher- 
ical geometry of the cluster is important for r > 1000 kpc. 
However, the model should represent the data fairly well 
within the range 500 < r < 1000. It is apparent from 
this figure that each of the three types of transport can 
produce the required fall off with distance, but that the 
associated predicted softening in the spectral index, which 
is also shown in the figure, is quite different in eac h case. 
The spectral index maps of Giovannini et al. ( 1993 ) , which 
were computed from the ratio of the fluxes at 326 and 
1380 MHz, show that within a few hundred kpc of the 
edge (at 500 kpc) the spectrum softens rapidly - by more 
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Fig. 4. The one-dimensional surface brightness profiles for the 
three types of transport a — 0.5, 1 and 1.5 fitted to observa- 
tions of the diffuse radio emission of the Coma cluster of galax- 



ies by Hanisch ( 198C ) as presented by Valtoja ( 1984 ) . Only the 
outer parts of the cluster are plotted. Also shown are the pre- 
dicted dependences of the spectral index on radius. The upper 
limit is a rough indi cation of the spectral softening observed 
by Giovannini et al. (1993) 



than 0.6 in 200 kpc. This is indicated by an upper limit 
in the figure. Of the curves presented in Fig. ^, only that 
corresponding to supra-diffusion produces a comparable 
effect. Thus, in the stationary state, and assuming in situ 
acceleration is negligible, diffusion is not capable of prop- 
agating the electrons outwards from the inner 500 kpc of 
the cluster. 

The conclusion that anomalous transport is necessary 
is based on assuming that both the transport coefhcients 
and the magnetic field (and also the type of transport 
as determined by a) are not functions of position, and 
further, that the transport is independent of particle en- 
ergy. Although it is not straightforward to include such 
effects in a transport theory, they do not provide a simple 
explanation of the observations shown in Fig. ^ in terms 
of standard diffusive transport. For example, a magnetic 
field which decreases with distance from the core, could 
explain the fall-off in the intensity, but could not account 
for the observed softening of the spectrum, since energy 
losses would then be less important at larger radius than 
assumed in the figure. Similarly, if the spatial diffusion co- 
efhcient increases in proportion to particle energy - as ex- 
pected in the case of gyro-Bohm diffusion - then particles 
responsible for higher frequency radiation are more mobile 
than those radiating at lower frequency. In this case, the 
intensity should fall off more slowly at higher frequency, 
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which would tend to produce a harder spectrum at larger 
radius, contrary to the observations. 

5. Conclusions 

In this paper we have presented a simple one-parameter 
method of modelling the effects of anomalous transport 
on energetic electrons. Anomalous transport is likely to 
occur wherever magnetised electrons (i.e., electrons which 
are tied to field lines) move in a magnetic field which has 
a stochastic component. In astrophysics, this is the rule 
rather than the exception and anomalous transport can 
be expected, for example, in the interstellar medium, as 
well as in the intra-cluster medium of clusters of galax- 
ies. The primary observational diagnostic is the intensity 
and spectrum of the synchrotron radiation emitted by the 
transported particles. For the simplest case in which the 
magnetic field is of constant magnitude and of random 
orientation, we have presented general expressions for the 
surface brightness as a function from position of injection, 
and also for the spatial variation of the spectral index. 

Applying these to the diffuse emission observed from 
the outer parts of the Coma cluster, we note that it is 
not possible to distinguish between the various forms of 
transport merely from the profile of the surface brightness. 
However, the expected spatial dependence of the spectral 
index is very sensitive to the type of transport. Assuming 
that the effects of particle acceleration are negligible in 
the outer parts of the cluster, and that the electron distri- 
bution has achieved a steady state, we find that standard 
diffusive transport cannot produce the observed rapid soft- 
ening of the spectrum with radius. Under these assump- 
tions, the type of transport indicated is supra-diffusion, in 
which particles move almost ballistically in a field config- 
uration which has an ordered radial component. 

The computations we have presented contain several 
major simplifications. In addition to the assumption of 
constant, randomly orientated magnetic field, and the sim- 
ple planar or spherical geometry, we have assumed that 
the parameters governing the transport are independent 
of the particle's energy. In reality, the type of transport 
itself (i.e., the value of a) will change depending on the en- 
ergy range and timescales considered. Thus, at very large 
times (which may exceed the synchrotron lifetime) , a par- 
ticle can be expected to decorrelate from the magnetic 
field and perform diffusion (e.g., Duffy et al., |1995D . We 
do not model the situation in which a significant change 
in a occurs within the synchrotron lifetime of an electron. 
Finally, in order to apply such models to well-observed ob- 
jects such as spiral galaxies, it will be necessary to include 
additional effects such as that of a galactic wind, as well 
as bremsstrahlung and ionisation loss processes. However, 
these processes will not change our basic conclusion that 
it is the spatial dependence of the synchrotron spectral 
index which provides the most sensitive measure of the 
transport properties of the emitting electrons. 
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A. Derivation of the propagators 

The propagators are evaluated by inserting the expansions 
(13) into Eq. (^0|), inverting the Fourier transformation 
analytically (since we fix /3 = 2) and using the method of 
steepest descents to find the asymptotic dependence of the 
S(^) at large ^. For the one-dimensional standard 
integral leads to 



P(x, s) = - exp 
s 



a/2 



(Al) 



with ^ = \/2{x/a){To/t)°'^^, from which, by the method 
of steepest descents (e.g., Mathews & Walker 1970| ) one 
finds 



cr v/27r(2 



exp 



2 — a \ I at, ^ 
Y 



-(1-q)/(2-q) 



(A2) 



In the case of diffusion (a = 1), this result is correctly 
normalised, i.e., J_^dxP{x,t) = 1. However, in gen- 
eral the normalisation is lost in this approach (cf. Rax & 



White 1992). Correcting for this, we find for the constants 



used in Eq. (17) the expressions 

'2-a\ /a\2/(2-a) 
.2. 



A 



1 



d{2-a) /a\"/(2-") 



: — aj /a\ 
V "2/ 



d/2 



(A3) 
(A4) 



B. The synchrotron brightness distribution 



The expression for the density (26) can be directly inte- 
grated numerically. Alternatively, it is easily computed by 
expanding the factor {1 — gjt')P~'^ and integrating over t', 
which yields 



n(a;,7) 



Qo7 



where 



i?(7) 



2xs ( a 



Vd\2-a 



and 



h„{n) = 1 , hk{n) = ■^^n(n - l)...{p - k + 1) 
fc! 
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Note that for p integer and larger than 2, the series re- 
duces to p — 2 terms. 

This expression recovers, for 3-dimensional diffusion 
(a = 1, d = 3), the form ula derived by Wilson ( 1975 ), and 
used by Valtaoja ( |1984D 



= 87r3/2jj3/2 



3 57r^ 



' AD 



2' AD 



(B2) 



where U is the confluent hypergeome tric fu nction of the 
second kind (Abramowitz and Stegun 1972 ). This can be 
seen either directly from Eq. (26) or, in the special case 
p = 2, by noting that 
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